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SUMMARY 

Presented is a description of a technique for the optimization of airfoil 
pressure distributions using an interactive inverse boundary-layer program. 

This program allows the user to determine quickly a near-optimum subsonic 
pressure distribution which meets hi? requirements for lift, drag, and pitching 
moment at the desired flow conditions. The method employs an Inverse turbulent 
boundary- layer scheme for definition of the turbulent recovery portion of the 
pressure distribution. Two levels of pressure-distribution architecture are 
used — a simple roof top for preliminary studies and a more complex four-region 
architecture for a more refined design. A technique is employed to avoid the 
specification of press«:'*e distributions which result in unrealistic airfoils, 
that is, those with negative thickness. The program allows rapid evaluation 
of a designed pressure distribution off-design in Reynolds number, transition 
location, and angle of attack, and will compute an airfoil contour for the 
designed pressure distribution using linear theory. 

Comparison of pressure distributions and corresponding airfoil geometries 
resulting from different specifications of recovery-region bound ary- layer form 
parameter will be presented. 


INTRODUCTION 

AirfbH design has traditionally been a trial-and-error process. Before 
large computers were available the airfoil designer had little option other 
than to make small changes to existing sections, guided by linear theory, in 
the hope that he would improve the section's performance in the area he wanted. 
Or, if he could afford the time and expense, the designer could undertake the 
testing of a large number of sections to try to find a good one for his appli- 
cation. Understandably, airfoil progress was a rather hit-or-miss proposition. 
Certainly, airfoil designers were guided by the physics of the flow around their 
sections but the mathematics required to describe the flow field adequately was 
too complex to solve in closed form and without computers, and it was too time 
consuming to solve numerically. 

When relatively cheap computing equipment became available, methods for 
solving the Inviscld flow fiel.d about an arbitrary airfoil evolved, and, soon 
after, methods for computing boundary-layer characteristics were developed and 
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Cftuiriii These were powerful tools for the airfoil designer. He 

It A * variations in his design variables and could be 

af had been previously tested. However, as good 

variabLrinhraJrfon"r\ ’'^***^^ nonlinear way the different 

and-trv wn? ‘ ^ P^^lem interact there is no guarantee that cut- 

to de^Ln i K? rr^ l««provement in performance; it is all too easy 

it i> JoLtir T t ? o“« aspect of performance only to find 

it is totally unacceptable in another. 

8et around the inherent problems of cut and try. the Inverse of the 
outiM^thrnlS* problem needed to ba colved more rigorously. Instead of com- 
the characteristics of a given section, one needed to specify 

oJobl!«^?rt"r compute the section. At present, this inverse airfoil 
problem is being attacked in two parts, as was airfoil analysis, that is the 
averse inviscid and the inverse boundary-layer problems. * 

and rli® fr®"® ^“^^scid problem, that is. specifying a pressure distribution 

by se^eJal XSleir been well developed recently as evidenced 

tL papers on the subject given in this volume. However, while 

iell in^hanr®i-^°*^ performing the inverse inviscid computations are seemingly 
off^ the use of these techniques for airfoil design is not so well 

tribuM^n fo^®" although there is a real, unique inviscid pressure dls- 

TtlTrlTJT ® ®"«^® the reverse is not true in 

J easily specify a pressure distribution which results 

iaHd negative thickness. Although such an airfoil is a perfectly 

liii? folution to the inverse inviscid p-oblem in a mathematical sense, one 
defiLd^''® ® ^ “i? trouble building one. The fact is that most arbitrarily 
ined pressure distributions result in unrealistic airfoil contours if present 

iTellTiclTilTof hlindly. Future programs which allow 

j f ^ pressure distribution and thickness by using welghtine 
and advanced solution techniques (such as least-squares or Newton's methodl^to 
resolve the over-specified nature of the problem Ly aid thHes^gn. 

iJ tL’^Ihoir® ?®'"® ^®^ in general use. So. one must be very careful 

are L iL^^o’^Id^n^ge!^^ '' 

How should the pressure distribution be defined? With a little thought 
d K®® numbers the performance of an airfoil is either 

(or buildable) thickness form. In order to proceed with a rational design pro- 

met^r® ® ® technique that, given boundary- laye r paJa- 

nrnd ’ compute a pressure distribution and, hopefully, is constrainL to 

produce only pressure distributions that are realistic. 

of thl^TirTT A tool in hand the designer might take advantage 

^ design process shown in figure 1. The most notable aspect of 

s process is that it proceeds from performance requirements to initial con- 

trlhuH il! *^*^® ^”''®’^®® ®°*^® " to compute a desirable pressure dis- 

trtbutlon and then to compute an initial airfoil contour. The last block. 

the detailed analysis and refinement stage, is the last remnant of trial and 
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distribution* alrfou from a pressure 
aerodynamics) the primary task of*^h a problem in mathematics*^ not 

'fr/“*ta'ier of thL rixn'jfr tor the job >t 

believe will greetly aid the designer ij JhlJ jjjt ‘ “'hnldue «hlch 


The synbols nsed herein are defined In an appendix. 


TECHNIQUt: 


The Inverse Boundary-Layer Equations 


equatinn^nhinn nnnnnnr® ® the boundary-Uyer 

trlbutlon Is 

^ pressure dls- 

used the boundary-layer momentum eqSatloi cll <lifcu88 in this paper, I have 
meter variation, and the Ludwis-TllW ’ ® equation for form para- 

(ref. 1). These equations^! beenT^rJ^^" shear™?nnn 

technique is shown as follows: ^ «mecer H ■ 6*/6. The solution 


saving Gamer's e, nation for the veXotlt, derivative glvea 


d(s/c) 


- dH/d(x/c) 


(u/u 


0.0135(H-1.4)(u/u ) 


(0/c) 


From the momentum equation solve for 


ilQ/c.) . Io_ + 2)(e/c) d(u/u )/d(s/c) 

d(s/c) ^ 2 “ — rr: — 
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Is given by Ludwlg-Tillman: 


The well sheer stress coefficient 


Tq/PU^ 


- 0.123(R-)“ X 

pu^ ® 


( 3 ) 


The momentum 
equations (1) 
dH/d(x/c)): 


velocity ratio are given by integrating 
and (2) numerically with a known H(x/c) (and thus a known 


0/c - 


<i(9/c) 

d(s/c) 


As/c 


(4) 


d(u/u^) 


(5) 


integration step in arc length. Equations (1) to (5) 
should be relaxed at each step for average values of u/u and 6/c. 


Pressure-Distribution Synthesis 

turbulent boundary-layer technique and an appropriate 

butloron ^ ” architecture one may quickly design a pressure distri- 

bution on one surface of an airfoil. To explain the way the pressure dlstri- 

architectures Lai^le - 

region startiL ft th* characterized by an acceleration 

® at the leading edge and terminating at an input fair point. 
Constant pressure is assumed from the fair point to the input beginning of 

desJrfd ti^t? f®cilltates pressure recovery from the roof-top pressure to the 
found bv thd» ”*"®*^*® P*^®®®'*’^®' The roof-top pressure is not input and Is 
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Iterative Procedure for Determining Minimum Cp 


Given a fair point, a recovery point, a trailing-edge pressure, and a guess 


at roof-top pressure C 


Pmin* 


the pressure distribution is assembled up to the 


recoveiy point. This pressure distribution is then analyzed to provide the 
starting values cf H and 0/c for the inverse turbulent boundary-layer module. 
Having these values, the recovery pressure distribution is computed from 
at the recovery point to the trailing edge by using the desired variation of 
form parameter in the recovery region. If the computed trailing-edge pressure 


is not that desired, 


is incremented, the pressure distribution up to the 


recovery point is reassemfe?ed, and the process is repeated to convergenc* 2 . 


The Lower Surface Pressure Distribution 

Sine 2 the total pressure distribution must represent a realistic airfoil 
and the upper surface pressure distribution is defined by boundary-layer con- 
siderations alone, the lower surface pressure distribution must be defined by 
the thickness requirements specified by the designer. In the present program, 
a standard th-f"*!:ness form is used, the NACA OOXX, which is scaled to the 
designer's desired maximum thickness. So with a single input t/cmax» the 
designer ^-^ts a lower surface pressure distribution which will result in an 
airfoil which has the upper surface shape required to give him his designed upper 
surface pressure distribution and a lower surface contour which results from an 
NACA OOXX thickness form cf the desired maximum t/c. Linear airfoil theory is 
used to accomplish this and I will not discuss it further in this paper. 


The Four-Region Architecture 

Although the roof-top architecture is effective and simple to use for the 
preliminary stages of pressure-distribution optimization, it allows very little 
flexibility for controlling the laminar and transitional portion of the boundary 
layer. As will be shown later, the transition-point position is often of first- 
order importance to the airfoil-design problem. To allow for more precise 
control of the laminar boundary layer, the four-region pressure distribution 
was devised. 

This architecture is shown in figure 3 and is characterized by an accelera- 
tion region (1), a region of constant pressure gradient (II), a laminar stress- 
ing region (III), and finally the turbulent recovery region (IV). The fair 
point is input as before; and, in addition, the desired value of pressure gradient 
is given for region II, and a desired point of initiation of laminar stressing 
and a desired value of laminar form parameter are given for region III. A 
simple, inverse, laminar boundary routine is used to produce the pressure dis- 
tribution required to produce the desired laminar forio parameter. Region III 
may be used for either stressing the laminar boundary layer to achieve rapid 
transition without separation or may be used to avoid transition - depending, of 
course, on the value of H specified. The equations used in this region are 
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u/u^ = C(s/c) 


/l - H/2.55^ 
\ 0.94 / 


C 


(u/u ) / (s/c) 

00 Q 


( l - H/2,55 \ 
\ 0.94 / 


where ( )^ refers to the values at the beginning of region III. 

In the derivation of these equations some liberty was taken in dropping 
terms for simplicity; however, practice has shown this relationship is remark- 
ably accurate. The symbol Cp . is defined as the pressure at the beginning 
of region III. Region IV, the turbulent recovery, is defined by the inverse 
turbulent boundary- layer technique as described earlier. 


Optimization Using an Interactive Program 

The pressure-distribution design technique I have described has been im- 
plemented in an interactive optimization program. That is, the program user is 
operating the computer program in a conversational mode where input is requested 
by typed messages and given to the program, real time, by responding with typed 
input at a terminal. Output is printed or plotted immediately at the terminal 
at the request of the user. 

Why an Interactive program? Pressure-distribution design is a highly over- 
defined optimization problem. Constraints and secondary requirements are tough 
to quantify and the weight given them often reflects the judgement of the 
designer alone. In deciding how to implement this technique I adopted the 
philosophy that it is desirable to leave the judgement in the hands of the user 
and that a computerized pressure ^distribut ion design tool should eliminate the 
computation in the design process and amplify the user’s judgement by quickly 
and clearly illuminating the important relationships between the variables and 
constraints of the design problem at hand. 

An interactive program is a natural outcome of this philosophy as it allows 
an adaptive flexibility in the optimization process impossible using batch-type 
computing. Furthermore, and not a bit less important, we have found that being 
a part of the optimization process is very educational for the user, especially 
when he is designing in a flow regime where he has little experimental experience. 
Lastly, it is one heck of a lot more fun to compute this way - as anyone who has 
tried it will agree. While this may seem a trivial concern - the fact that the 
computation is now pleasant, rather than a chore, often means one will tend to 
stick out a really tough problem longer. 

The options available ia this optimization program are illustrated in 
figure 4. I will just summarize them here. 
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design is Initiated by giving the program the flow conditions and 

^hef form-parameter variation is 

then chosen - this program allows a constant form parameter or linear, quad- 

arM^rarv variation to be specified with very little input. An 

arbitrary H variation may be used if a file of th- desired values of H vs. 

T generated. This complete, the program enters the 

?>.r J ^ dis- 

tribution. This whole process takes about 30 seconds. At this point Me user 

Stlsfie^with he any other of the design paramecers. or if he is 

satisfied with Lhe design-point pressure distribution he may analyse the pressure 

istribution at an off-design Reynolds number or trip location or ne may analyse 
the pressure distribution off-design in lift (accomplished by adding an angU- 
attack velocity distribution to the design distribution using linear theory), 
he results of any of these off-design analyses are printed and plotted 
Wdlately at the terminal. In the event that a desirable pressure distri- 
on is generated, c nrellmlnary airfoil contour may be obtained simply by 
request. This is accomplished by a simple linear inverse program due to^ 

airfoils produced should be considered as starting 
points for more accurate design methods. They are, however, excellent bench 
mrks to check the resulting airfoil's fidelity to structural restraints (such 
as a restraint on maximum camber or compound surfaces) • 

‘^J® P^®^ ‘^his approach, consider that to define a pressure 
should be checked at several Reynolds numbers and off-design 

P>^®duce a preliminary airfoil takes about 7 minutes on 
?. .Y - or about 8 minutes on a CDC Cyber 175. The reason the more power- 

^nn *^*^® ‘^ata-transmlssion rate to our terminal is 

r *'*'® time is used during 

transmission and user head scratching in deciding what to do next. ® 


Some Illustrative Results 

® comparison of pressure distributions and their corres- 
5 geometries resulting from different specifications of recovery- 
region boundary-layer form parameter (H) . The 'C, 'L', and '£' deiote con- 
stant H, linearly varying and exponentially varying H, respectively. Each 

?hrimnL^ designed to a lift coefficient of approximately 1.6. 

he impacc of the variation of the form parameter on both airfoil performance 

99 «; f striking. The L/D (e*) ranges from 172 for constant H (2.06) to 

225 for the exponentially varying H (1.45 at x/c = 0.3 to 2.0 at x/c = 1.0) 
Pltc.ilng moment varies from 0.057 to -0.186 for the same two examples. The 
airfoil contours, I think, speak for themselves. Again, the values of H for 
each type of variation were picked so that the pressure distributions all pro- 
duced a design-point lift coefficient of about 1.6. 

In the designs of figure 5, natural transition occurred near the fair 
point (x/c = 0.1) due to very high design-point Reynolds number (Rn - 30 x io6). 

J illustrate the powerful effect that transi- 

tion has on the design of a pressure distribution. In figure 6 a design was again 
undertaken at a Reynolds number of 30 x iq 6 with a roof-top archltecttife and 
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natural transition. The result was a pressure distribution of modest perfor- 
mance with transition occurring very near the fair point. 

Using the four-region architecture, an attempt was made to maintain laminar 
flow as far aft on both surfaces as possible. The result was a pressure dis- 
tribution that produces significantly better performance. The L/D at the 
design point went from 215 for the roof top to 317 for the four-region Press 
distribution. Design was Increased from 1.42 to 1.68 for 

fication in form parameter, whereas the pitching-moment coefficlen 
relatively unchanged. Although the final pressure distribution 
shows remarkable improvement and impressive performance, work still needs to 
done to achieve acceptable off-design performance as one might expect considering 
the section's rather sharp leading edge. 

Figure 7 displays the effect of tripping the laminar boundary layer on 
design lift and L/D at a Reynolds number of 4 x 10 for several recovery oca- 
tlonL Laminar flow is extensive when no trip is specified due to the lower 
Reynolds number, even though the roof-top architecture is used throughout, so 
one would expect tripping to have a significant impact. The curves 
does. The peak design is reduced from 2.25 to 1.3 oy tripping and the L/D 
is reduced from 177 to 77. The explanation for this sharp reduction in 
formance is that early transition produces a much thicker boundary layer at the 
recovery point due to the more rapid growth of a turbulent boundary-layer rela- 
tive to the laminar layer. This thicker boundary layer is able to withstand 
less adverse pressure gradient and so must recover from a lower velocity. Thus, 
there is a loss in lift and L/D. 

Figure 8^ summarizes the performance of some of the pressure distributions 
we have designed to date. With the exception of points 3 and 4 all were low C^, 
designs with heavy restrictions placed on the maximum camber and all the 

simple roof-top type architecture with tripping to produce a turbulent bounda y 
layer at the recovery point where necessary. Because of this, one "oj 

interpret the results as an attempt to plumb the limits of L/D at each Reynolds 
number but rather consider these results only as an example of the type o 
information obtainable from this program with a relatively modest amount of time 
and expense. Points 3 and 4 were attempts at optimizing L/D at their design 
Reynolds numbers, and although they do perform credibly, more effort is required 
to refine the designs and guarantee good off-design performance. 


CONCLUSIONS 


We need more experience with this program to say with any certainty 
whether it can be used to determine an optimum pressure distribution for t e 
flow conditions where the theories are valid. However, the following aspects 
of this technique are already apparent: 


^Dr. John McMasters assisted in preparing this paper and allowed me to use 
the data in figures 5 and 8 which he generated for a Boeing internal study of 
the L/D potential of moderately cambered airfoils. 
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1. The Inverse boundary- layer technique is a powerful and efficient device 
for defining a pressure distribution. 

2. The pressure— distribution architecture need not be excessively complex 
to control the laminar and transitional regions of the pressure distribution. 

3. The simple Inverse laminar equation is an effective tool for stressing 
the laminar boundary to transition without causing separation. 

4. Allowing a scaled thickness form to define the lower surface pressure 
distribution does constrain the program to produce pressure distributions which 
result in realistic sections. However, more flexibility is needed in this area 
so that aerodynamic considerations may also play a part in definition of the 
lower surface distribution. This may take form as optional thickness forms or 
arbitrary thickness forms constrained only in maximum thickness and to positive 
thickness everywhere. 

5. The interactive type of execution seems to have value beyond the obvious 
advantages of speed and immediate response in the display of Intermediate re- 
sults, that is, it is very educational and pleasant to use. 
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APPENDIX 


SYMBOLS 


H 

U 

U 

00 

s/c 

x/c 

Rn 

"e 

0/c 


6*/c 


P 

C 

C 

c 


£ 

d 


m 

C 


P 


t/c 




c 


boundary-layer form parameter, 6*/0 
velocity at outer edge of boundary layer 
free-stream velocity 

airfoil-surface arc length normalized vfith chord 
distance along airfoil chord, normalized with chord 
Reynolds number based on airfoil chord 

Reynolds number based on local value of momentum thickness 6 and 
velocity U 

boundary-layer momentum thickness (normalized) , 



^ (1 - ^) 


dy/c 


boundary- layer displacement thickness. 



(1 - -^) dy/r 

oo 


density 

section lift coefficient 

section drag coefficient 

section pitching-moment coefficient 

pressure coefficient, C = 1 - 

P 

airfoil thickness (normalized) 
design-point section L/D 
maximum section L/D 



I 
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Figure 1.- Airfoil design process. 


End of Initial Acceleration (Fair Point) 


Beginnintj of the Turbulent 
Recovery (Recovery Point) 
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Form Parameter Va tion 



x/c 


Trailing Edge Cp 




Lower Surface Cp Determined 
by Desired t/c Max 


Figure 2.- The roof-top architecture. 
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Modify Cp Distribution for Changed Lift 
Coefficient by Adding Angle of Attack Distribution 


Figure 4.- Block diagram of interactive 


pressure-distribution optimizer. 
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BHL 120-3OO-iO/?B IH Variation) 


H 

Variation 

ci 

Cm 

€* 

Crux 

ClLiu 

E 1.45/2.4 

1.592 

-0.186 

226 

24C 

2.3 

L 1.45/2.0 

1,563 

-0.058 

216 

216 

2.05 

C 2.06 

1.603 

0.057 

172 

172 

1.97 


Figure 5.- Comparison of pressure distributions and corresponding airfoil 
geometries resulting from different specifications of recovery-region 
boundary-layer form parameter H. 


BHL 120- 300-40/25-El. 45/2. 4 

C* • 1.42 C* - 0.00661 Cd - 0.00520 
I a “U 

C* - -n.214 c* » 215 


BHL 120-300-40/38-El. 45/2. 4 HOD, 

C* * 1.68 cj - C.U353 - 0.0051 

C* • -0.237 C* - 317 “ 






Figure 6.- Examples of basic pressure-distribution architecture. 


396 





